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ABSTRACT 

Context. The measurement of the anisotropies in the cosmic infrared background (CIB) is a powerful means of studying the evolution 
of galaxies and large-scale structures. These anisotropies have been measured by a number of experiments, from the far-infrared to 
the millimeter. One of the main impediments to an accurate measurement on large scales (<1 degree) is the contamination of the 
foreground signal by Galactic dust emission. 

Aims. Our goal is to show that we can remove the Galactic cirrus contamination using Hi data, and thus accurately measure the 
clustering of starburst galaxies in the CIB. 

Methods. We use observations of the so-called extragalactic ELAIS Nl field at far-infrared (100 and 160 /um) and radio (21 cm) 
wavelengths. We compute the correlation between dust emission, traced by far-infrared observations, and Hi gas traced by 21 cm 
observations, and derive dust emissivities that enable us to subtract the cirrus emission from the far-infrared maps. We then derive the 
power spectrum of the cosmic infrared background anisotropies, as well as its mean level at 100 iim and 160 iim. 
Results. We compute dust emissivities for each of the Hi- velocity components (local, intermediate, and high velocity). Using 
IRIS/IRAS data at 100 //m, we demonstrate that we can use the measured emissivities to determine and remove the cirrus contri- 
bution to the power spectrum of the cosmic infrared background on large angular scales where the cirrus contribution dominates. We 
then apply this method to Spitzer/MIPS data for 160 iim. We measure correlated anisotropies at 160 yum, and for the first time at 
100 fim. We also combine the Hi data and Spitzer total power mode absolute measurements to determine the cosmic infrared back- 
ground mean level at 160 fim. We find Z?i 6 o = 0.77+0.04+0.12 MJy/sr, where the first error is statistical and the second one systematic. 
Combining this measurement with the B\ o/B\eo color of the correlated anisotropies, we also derive the cosmic infrared background 
mean at 100 jim, Z?ioo = 0.24+0.08+0.04 MJy/sr. This measurement is in line with values obtained with recent models of infrared 
galaxy evolution and Herschel/PACS data, but is much smaller than the previous DIRBE measurements. 

Conclusions. The use of high-angular resolution Hi data is mandatory to accurately differentiate the cirrus from the cosmic infrared 
background emission. The 100 fim IRAS map (and thus the map developed by Schlegel and collaborators) in such extragalactic fields 
is highly contaminated by the cosmic infrared background anisotropies and hence cannot be used as a Galactic cirrus tracer. 
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1. Introduction 

Starburst (SB) galaxies are known to have had an important role 
in galaxy formation and evolution throughout the whole history 
of the Universe. In the far-infrared (FIR) and sub-millimeter, 
observations are limited by extragalactic confusion: details on 
small spatial scales are lost in the noise because of the poor angu- 
lar resolution of the instruments. As a result, unresolved starburst 
galaxies form the cosmic infrared background (CIB) ( [Puget et aL] 



resolved -45% and 



1996; Fixs en et al.][T998| |Lagache et al. 1999), which peaks 



at around 200 fim. In the mid-infrared, a large fraction of the 
CIB has been resolved into individual sources: |Papovich et al.| 
(2004), for instance, resolved 70 % of the 24 /im background. In 
the FIR, before the advent of the Hersch el telescope , a sm aller 
fraction has been resolved: with Spitzer, [Dole et al.| ( 2004| ) re- 



solved 23% and 7% of the CIB at 70 and 160 yum, respectively. 
Frayer et al.] ( |2006b| ) managed to resolve 60% of the CIB at 70 
fim using a very deep but small field, which was thus limited 
by cosmic variance. More recently, Berta et al. ( 2010 ) integrated 
counts coming from Herschel/PACS data at 100 and 160 yum and 



wavelengths, |01iver et aL 



52% of the CIB, respectively. At longer 
" 20T0]) directly resolved 15%, 10%, 



and 6% of the CIB at 250, 350, and 500 //irrespectively, using 
Herschel/SPIRE data. Confusion can be circumvented by the use 
of statistica l methods. For instance, by stacking 24 yum sources, 
Dole ~et~aLl ([2006) were able to resolve a large fraction of the 
CIB at 70 fim and 160 yum and |Berta et al.| ( [20T0| ) increased their 
fractions from 45% to 50%, and from 52% to 75% at 100 fim 



and 16 /mi, respectively. Using a P(D) approach, Bert a et al. 
( |2011| ) were able to obtain still larger fractions of 65% and 89%, 
at 100 and 160 yum, respectively; P(D) derived counts of Glenn 
|etafl { |20l0) account for 64, 60, and 43% of the CIB at 250, 350, 
and 500 /mi, respectively. These results imply that the sources 
detected at 24 fim constitute the bulk of the CIB around its peak. 



Caputi et al. (2006) showed that galaxies that dominate the emis- 



sion at 24 fim become more and more luminous and massive as 
the redshift increases starting from luminous infrared galaxies 
(LIRGs) with 10 n L o < L IR < 10 12 L o at 0.8 < z < 1.2 with in- 
termediate mass, to ultra-luminous infrared galaxies (ULIRGs) 
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1O 14 L that dominate at z > 2 and have 



with 10 12 L o < Lir 
masses > 1O U M . 

The clustering of galaxies that make up the CIB can be char- 
acterized by its anisotropies. This clustering was first detected 



at 160 yum with Spitzer (Lagache et al. 2 007} |Grossan & Smoot 



2007 ), and then measure d at 250, 350, and 500 /im using BLAST 
data ( |Viero et aT1|2009) . All three of these sets of data enabled 
the detection of an excess of signal on intermediate spatial scales 
caused by the clustering of starburst galaxies which enabled 
them to derive the linear bias parameter b that relates the density 
fluctuations of luminous matter to those of dark matter (DM). 
Lagache et al.|(|2007) measured a value of b = 2.4 ± 0.2 while 
Viero et al.| (|2QQ9 ) obtained b = 3 ± 0.3. The difference may 



be due to selection effects. At longer wavelengths, higher red- 
shift SB galaxies a re probed ( Lagac he" et al . | ( [2505] ) , |Fernandez- 1 
Conde etaL] ( |2008| )) and at these higher redshifts, SB galaxies 



are a highly biased tracer of the underlying dark matter density 
field. They indeed formed in very massive DM halos early in the 
history of the Universe. Magliocchetti et al. (2008 ) derived the 
two-point correlation function of 24 fim selected sources divided 
into two redshift bins (0.6 < z < 1.2 and z > 1.6), finding that 
these SB galaxies are strongly clustered and embedded in DM 
halos of 1 O 12 - 8 M for the high z sample and 1O 1L8 M for 
the low z one. |Cooray et al.| ( |2010) computed the angular correla- 
tion function with Herschel/SPIRE data. They found that 250 jim 
sources are embedded in DM halos of ~ 10 12 M o at < z >~2.1, 
whereas bright 500 /jm sources reside in more massive halos, 
~ 10 13 M o , at < z >~2.6. The CIB anisotropy measurements in 



the F IR and submillimeter were followed by those of |Hall et aL 
( 2010) at 1.3 mm and 2 mm with the South Pole Telescope a nd 



by the Fowler & Atacama Cosmology Telescope Team|(|2010|) at 



1 .4 and 2 mm with the Atacama Cosmology Telescope (Dunkley 



& Atacama Cosmology Telescope Te am (2010)). More recently, 
Planck Collabora tion et al.| ( |2011b| ) derived CIB power spec- 



tra from 10' to 100' simultaneously at 3 50, 550, 850, and 1380 
/zm in six high-Galactic latitude fields. Ambl ard et al. ( 2011| ) 
extended the measurements to smaller angular scales, using 
Herschel/SPIRE observations of the Lockman-hole field at 250, 
350, and 500 /zm. These measurements allow us to start to refine 
the analysis of the clustering properties of galaxies responsible 
for the CIB, and its cosmic evolution to high redshift (z~3-4). 
The far-infrared and submillimeter emission of the Galactic cir- 
rus interferes with the detection and measurement of the CIB. 
This emission dominates the power spectrum of the anisotropies 
on large spatial scales, the exact scale depending on the wave- 
length and the selected field. To remove the cirrus contribution, 
data at other wavelengths are usuall y used. For instance, IRIS 
maps of reprocessed IRAS maps ( Miville-Deschenes & Lagache 
2005 ) at 100 jim can be used to determine the power spectrum of 
the cirrus on large scales. However, as shown in this paper, these 
data also contain CIB anisotropy and therefore a clustering sig- 
nal from SB galaxies. To constrain more accurately the contribu- 
tion of the Galactic cirrus, an external tracer is needed. The most 
effective one for dust emission in the diffuse sky is neutral hydro- 
gen. In this paper, we show that Hi data can be used to remove 
the cirrus contamination from the 100 //m and 160 //m maps in 
order to measure the CIB intensity and the power spectrum of the 
CIB anisotropy. This method based on template removal to sep- 
arate the cirrus and CIB components, was also successfully ap- 
plied to the -140 square degrees of the very diffuse high-latitude 
sky observed by both Planck/HFI (350 /mi to 3 mm) and the 
Green Bank Telescope (21 cm Hi) (Planck Collaborat ion et al. 
|201 lb) . 

The paper is organized as follows : we present the data in Sect. 2. 



MIPS destriper 
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non destriped 
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Fig. 1. The dotted black line shows the power spectrum of the 
non destriped map, the dashed red curve is the power spectrum 
of the map destriped with the "IRIS destriper", and the green line 
shows the power spectrum of the map destriped by "the MIPS 
destriper". The vertical dashed line shows our angular scale cut 
of < 0.8 arcmin -1 . 




Fig. 2. The final source- subtracted 160 /mi map centered on 
ELAIS Nl used to compute the power spectrum. Units are 
MJy/sr. 



In Sect. 3, we briefly describe the as trophy sical components and 
derive the instrumental noise component of the map power spec- 
trum. In Sect. 4, we compute the Galactic cirrus contamination. 
The cirrus map is then subtracted from the IR maps to obtain the 
CIB -dominated maps from which we then estimate the power 
spectrum at both 100 /mi and 160 /mi (in Sect. 5). In Sect. 6, we 
estimate the mean of the CIB at 160 //m using the total power 
mode of Spitzer, and at 100 /mi using the #ioo/#i60 color of the 
measured correlated CIB anisotropies. We then present our con- 
clusions in Sect. 7. 



2. Data 

We focus our analysis on the ELAIS Nl field (£,b)= 
(85.33°,44.28°), which is part of the Spitzer Wide- Area Infrared 
Extragalactic Legacy Survey (SWIRE). It covers about ten 
square degrees and was observed by the Multiband Imaging 
Photometer for Spitzer (MIPS) at 160 /mi, by the GBT at 21 
cm and by IRAS at 100 and 60 /mi. 



2 



Aurelie Penin et al.: An accurate measurement of CIB mean and anisotropics 





Fig. 3. Left: point spread function (PSF) of MIPS 160 yum com- 
puted using an oversampled map. Right: PSFs power spec- 
tra. The continuous line is our PSF and the dashed line is th e 
power spectrum from the STiny Tim Model (Krist et al. (2005 ), 
STinyTim, vl.3; Krist, 2002). 



2.1. MIPS 160 nm 

The MIPS observations at 160 /im were take n as part of the 
SWIRE legacy program ( [Lonsdale eraT][2003] ) in 2004 during 
two epochs separated by six months (late January and late July). 
The data were taken using the scan mode and the medium scan 
rate, in about three-degree long strips and an offset between the 
return and forward scans of 148". This observing strategy pro- 
duced a map with a median depth coverage of 8 (times 4 seconds 
per frame = 32 sec integration time). However, owing to the dead 
readout of the 160 fim array ( Stansberry et al. 2007) the cover- 
age was inhomogeneous, being as low as two or as high as ten 
in some regions of the map, usually along the scanning direction 
and overlapping independent regions. 

During the observations of the first epoch, the Spitzer 
Observatory went into stand-by mode (Jan 25th). This meant 
that some of the 160 /mi observations were affected by a slightly 
warmer telescope during the recovery phase, with a temperature 
of around 6 K rather than the standard ~ 5.6 K. We mitigated 
this effect in the data by applying a small offset determined from 
the nearby unaffected regions (an overlap correction). 
The fact that the 160 /im data were taken during two different 
epochs makes the processing and creation of the final mosaic 
relatively straightforward. We used the standard basic calibrated 
data (BCDs) from the Spitzer Science Center, and removed their 
prediction of the zodiacal light as a function of time a nd space 
from each BCD and carried out an overlap correction (Mizuno 
et al.|2008 ). The maps were created using the native pixel scale 
of 16 arcsec/pixel, which preserves the diffuse emission. These 
final mosaics show no evidence of any defects other than low 
intensity stripes left as an artifact of the scanning observation. 
These stripes introduce two peaks at k ~1 arcmin -1 into the 
power spectrum of the map, as shown in Fig. [T] The effect of 
these stripes was mitigated by applying a destriping filter using 
ridgelets (Ingalls et al. 201 1). Our tests of the photometry of the 
maps and their sources showed that the destriping method pre- 
serves the flux to within 5% of the original values. However, the 
power spectrum of the destriped map shows a little loss of power 
on small scales (k>0.2 arcmin -1 ), as can be seen by comparing 
the green line with to the black dotted line in Fig. [T] To inves- 
tigate whether any astrophysical information had been removed 
or only power contained in the stripes, we destriped the map us- 
ing the destriper algorithm developed for IRIS, which has been 
shown to prevent the remov al of astrophysical signal (Mi ville- 
Deschenes & Lagache 2005). 

The power spectrum of that map is shown by the red curve 
in Fig. [T] It is in good agreement with the green curve up to 
k ^0.8 arcmin -1 , suggesting that the "ridgelet destriping" did 
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Fig. 4. The complete data sequence of one photometric observa- 
tions of the ELAIS Nl field at 160 /mi, with the median values 
of each frame in surface brightness units (MJy/sr) as a function 
Data Collection Event (DCE) number (essentially time). The 
highest values correspond to the calibration stim flashes. 
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Fig. 5. A zoomed image of the first 20 DCEs of the Photometic 
observation, showing the small exponential decay after a stim 
flash lasting about 2-3 DCEs. 



not remove any astrophysical signal. At higher k, we found 
that the IRIS destriper fails to remove residual stripes as it 
was not developed for such high-angular resolution data. We 
therefore use the map destriped with ridgelets. To limit point 
spread function (PSF) uncertainties (see below) we consider 
only scales <0.8 arcmin -1 , hence disregard the residual stripe at 
k ~ 1 arcmin -1 . 

The identification of sources and their extraction from the 
Spitzer long- wavelength images (70 and 160yum) has been car- 
ried out on filtered mosaics using the standard methods applied 



to the 'deep fields' (se e e.g.jFrayer et aL] ( |2006a) , |Frayer et al.| 
(2006b), Fray er et al. ( 2009| i). Filtering has been performed to 
remove the extended emission, thus enhance the detectability 
of faint point sources. Since we wish to preserve the back- 
ground/foreground emission, we carried out the source removal 
on our original mosaics. The PSF is well-sampled in MIPS 
data: the full width halh maximum (FWHM) of the 160 jim 
PSF is 40" compared to the pixel size of 16". We thus selected 
Star finder, which uses a PSF fitting algorithm, to perform 
our source extraction (Diolatti et al. 2000). Star finder 
also has the advantage that it evaluates the background over 
the entire image, in the course of the iterative fitting of the 
individual sources: this leads to smaller residuals (Marleau et al.l 



2004), even compared to those o f the standard Spitzer sourc e 
extraction software M0PEX/APEX flMakovoz & Marleau|[2005| ). 
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Fig. 6. As Fig.|4| but for the total power mode observation. Once 
again the largest values correspond to the stim flashes, with the 
sky measurement being shown in-between them. 
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Fig. 7. A zoomed image of the first 20 DCEs of the total power 
mode observation, showing as well the exponential decay, the 
sky measurement (DCE=8), small decay, and the next calibra- 
tion stim. 



We detected sources down to S i6o = 25 + 5 mJy and removed 
them from the map. Even though the source catalog is highly 
incomplete at low fluxes, removing sources down to very faint 
fluxes allows us to lower the Poisson part of the power spectrum 
and accurately measure the CIB correlated part. Fig. [2] shows 
the source- subtracted map that we use in our study. 

The effective PSF of the observed data is quite sensitive to 
the observed field. We therefore had to determine it directly in 
the map rather than use a prior estimate. We computed the PSF 
by stacking the brightest sources of an oversampled map (7.2 
arcsec/pix). The oversampling was necessary to obtain an accu- 
rate profile of the PSF. We extracted sources at 5<x and retained 
only those with S i6o > 600 mJy. The left panel of Fig. [^dis- 
plays the extracted PSF and the right panel i ts power spectrum 
compared to that of the STiny Tim Modef] Kr ist et al.| ( [2005] ), 
STinyTim, vl.3; Krist, 2002). The discrepancy between the 
power spectra shows the need to determine the PSF in our data. 
We estimated the errors in the PSF, using different oversam- 
plings, bright source flux cuts, and other extragalactic fields ob- 
served with the same scanning strategy. We found that our mea- 
surement is highly reproducible up to k ^0.8 arcmin -1 . At higher 
k, we have differences that can be as large as 35% which pre- 



1 http ://s sc . spitzer.caltech.edu/dataanalysistools/ 
tools/contributed/general/stinytim/ 



vented us from accurately measuring the CIB power spectrum. 
We will thus analyze the power spectra only for k <0.8 arcmin -1 . 

2.2. MIPS 160 jim total power modes 

In our analysis, we used of total power mode (TPM) MIPS ob- 
servations. A TPM observation at 70 yum and 160 yum is merely 
another example of a photometry mode observation in terms of 
its processing from raw data to calibrated images (basic cali- 
brated data or bed). As in the case of photometry, the calibration 
was directly tied to the stim flashes that are measured approxi- 
mately every two minutes. In the standard photometric observa- 
tions, the scan mirror places the source/sky on different positions 
of the array (s), taking an image of a frame at each correspond- 
ing position. Each frame is calibrated using the stim flashes, and 
goes through several steps as described in Gor don et al.|p005] ), 
described by equations (4) to (16), with two steps being funda- 
mental and applied to all the source frames, a dark subtraction 
and an illumination correction, i.e. division by a sophisticated 
flat. The total power mode differs in one fundamental way from 
the photometry mode, in that instead of using the scan mirror to 
move only the source around the array, it performs an interme- 
diate step of placing the mirror on an internal dark, an 'absolute 
reference' frame, between the source observations. This internal 
reference dark is inside the cold instrument, i.e. a temperature of 
less than 1.5 K (recall that the Spitzer telescope cryogenic min- 
imum temperature was 5.6 K). Laboratory and flight tests have 
shown that the attenuation is at least a factor ten when the scan 
mirror was placed at the 160 yum dark position. This 'reference' 
frame is observed in the same way as the sky or source frames 
i.e. is calibrated using stim flashes and undergoes a dark subtrac- 
tion and illumination correction. As pointed out by Gordon et al. 
(2005) in their Sect. 5.2.1, one problem when flashing a stim on 
the 70 /im and 160 yum Germanium arrays is that of a memory ef- 
fect ('stim flash latent') soon after the flash. This latency decays 
exponentially in a way that depends on the background light, so 
it can last 5-20 s and has a peak amplitude of less than 3% for 70 
jim and 7% for 160 yum. Experiments have shown that for signals 
of a few MJy/sr, similar to those in this study, the latencies are al- 
most completely reduced. We note that TPM observations were 
meant to be performed on diffuse signals and relatively dark re- 
gions of the sky, hence the design of the observing mode was 
done by taking into account this memory or latency effect. A 
standard 160 yum TPM observation places the array between the 
sky and the 'internal reference' (aka 'internal dark'), but after 
each observation the memory effect is allowed to decay by taking 
a series of seven frames, and using the frames unaffected by the 
latency of the final measurement. Figures [4] through [7] show the 
time history, as a function of the frame or Data Collection Event 
(DCE) number, of the median value for each BCD (i.e. the me- 
dian over the 2x 16 160 yum array) from one of photometry and 
TPM observations carried out at 160 yum for the ELAIS Nl field. 
These observations were done in a sequence, i.e one after the 
other. Fig.|4]and|5] show the photometric observation; the high- 
est values -50 MJy/sr corresponds to the stim flashes and the 
rest the sky measurement. The zoomed image shows the first 13 
DCEs in-between the two stims (DCE=1 and DCE=15) approx- 
imately two minutes apart. We note how pairs of frames after 
the stims are affected by the latency memory effect, and that the 
overall level is ~ 4.5 MJy/sr. Fig |6]and|7] show the TPM observa- 
tion, with once again the median values of the individual BCDs 
as a function of time (DCE number). The strongest signal is that 
of the stims and the weakest signal that of the source/sky itself. 
The zoomed image shows the first 20 DCEs, where DEC=0 cor- 
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Fig. 8. ELAIS Nl/GBT column densities N HI , in unit of 10 20 atoms/cm 2 . From left to right: the column density of the local compo- 
nent with -14 km/s < V LSR < 43 km/s, the IVC with -79 km/s < V LSR <-14 km/s and the HVC with -163 km/s < V LSR <-79 km/s. 
The white contour in each image shows the field covered by Spitzer (ELAIS Nl/MIPS field). 



responds to the sky measurement, DCE= 1 the stim, DCE= 2 to 
7 to observations in the "internal dark reference" and DCE=8 is 
the sky measurement followed by another six DCEs and the next 
calibration stim (DCE=15). We note that the level of the DCE=7 
is essentially zero, because in the data reduction pipeline this is 
the 'reference' measurement that is subtracted from each BCD. 
We also note that the sky itself (DCE=8) is at ~ 3 MJy/sr, which 
is the final TPM measurement. The comparison of the TPM with 
the photometric observation shows that the light contribution of 
the telescope at 160 yum is ~ 1.3 MJy/sr. The level of spurious 
emission due to the telescope background emission at 160 fim is 
1.0+0.2 MJy/sr. 

The absolute calibration of 160 yum TPM relies on the standard 
160 /mi calibration, which is based on asteroids, that is tied it- 
self to the 24 and 70 fim MIPS absolute calibration to be inter- 
nally consistent ( Stansberr y et al.|2007| ). Finally, we stress that 
a single DCE at 160 fim does not cover the 160 yum beam (40"), 
hence the TPM mode was designed to move the scan mirror to 
cover the beam in one cycle. A standard TPM 160 yum observa- 
tion contains four cycles. The final product for a single 160 yum 
TPM observation is a 5'x5' small mosaic, and for our mea- 
surements we calculated the mean over such an image. The two 
ELAIS Nl TPM observations at 160 jim took 2386 s with 88 
s on source each, while the two photometric observations using 
the enhanced mode took 671 sees with 54 sees on each source. 
These numbers illustrate the efficiency of the two modes; for ev- 
ery second on source at TPM, about seven seconds woth of data 
are used for calibration and latency decay. 



2.3. Hi data 

We used the Hi 25 deg 2 data cube (x, y, velocities) centered on 
ELAIS Nl. These data were obtained in 2006 and 2010 with 
the 100-meter Green Bank Telescope (GBT). Spectra were mea- 
sured over a 5° x 5° area centered on (£,b) = (85.5°, +44.3°) 
every 3.5' in both coordinates. Data were taken by in-band 
frequency switching yielding spectra with a velocity coverage 
-450 < Vlsr ^ +355 km s" 1 and a velocity resolution of 
0.80 km.s -1 . Spectra were calibrated, corrected for stray radi- 
ation, and placed on a brightness temperature (7&) scale as de- 
scribed in |Blagrave et al.| ( |2010| ) and |Boothroyd et al.| ( [20TTb| ). A 
third-order polynomial was fit to the emission-free regions of the 



spectra to remove any residual instrumental baseline. The final 
data cube has a root mean square (rms) noise in a single channel 
of 0.12 K of Tb, and an effective angular resolution of 9.4' x 9.1' 
in I and b, respectively. 

We distinguish three velocity components in the Hi gas data: 
the local, intermediate (IVC), and high velocity cloud (HVC). 
These are shown in Fig. [8] The HVC is centered around -115 
km/s and the IVC around -23 km/s, as illustrated in Fig. [9] This 
figure shows velocity spectra along three lines of sight, each of 
which is dominated by one component. The IVC and the HVC 
are clearly seen in the middle and bottom panels. 



2.4. IRIS/IRAS data 

We used IRIS (re-processed IRAS data) maps at 60 and 100 yum 
to measure the emissivities of the dust correlated to the Hi com- 
ponents and derive the CIB power spectrum at 100 yum. This 
new generation of IRAS images was processed using a more re- 
liable zodiacal light subtraction, from a calibration and z ero level 
compatible with DIRBE and a more reliable destriping (Miville- 
|Deschenes & Lagache||2005| ). At 100 /mi, the IRIS product 
also represents a significant improvement on the Schlegel et al. 
( |1998| ) maps. IRIS keeps the full ISSA resolution, it includes 
well-calibrated point sources, and the diffuse emission calibra- 
tion on scales smaller than one degree was corrected for the 
variation in the IRAS detector responsivity with scale and bright- 
ness. 

Using IRAS, two full-sky maps (HCON-1 and HCON-2 for 
hours confirmation) were performed and a last one that covers 
75 % of the sky (HCON-3). The three of them were processed in 
the same way, including deglitching, checking of the zero-level 
stability, visual examination for remaining glitches and artifacts, 
zodiacal light removal, and gain calibration. The three HCONs 
were then coadded using sky coverage maps to produce the av- 
erage map (HCON-0). We return to these HCONs later to deter- 
mine the power spectrum of the instrument noise. 
The IRIS PSF is assumed to be Gaussian following |Miville- 
peschenes et al.| ( |2002| ) 



P(k) = e X p\- — 



(1) 
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MIPS 160 jim IRIS 100 fim IRIS 60 /mi GBT 21 cm 



FWHM (arcmin) 



0.64 



4.3+0.2 



4.0 +0.2 



9.1 



Table 1. Effective measured PSF full width half at maximum (from 
and|Boothroyd et aL|(|2011a|) at 160 //m, 100 /mi and 21 cm, respectively) 



Stansberry et al. (2007), Miville-Deschenes & Lagache 



(2005) 
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Fig. 9. Hi spectra of three different lines of sight, illustrating the 
three velocity components (local, IVC, and HVC from the top to 
the bottom panels, respectively). 




?.3 



Fig. 10. Source-subtracted IRIS 100 yum map centered on 
ELAIS Nl. The area is ~ 25 square degrees. Units are MJy/sr. 
The white contour shows the MIPS field. 



3. Power spectrum and error bars 

There are several contributions to the power spectrum measured 
in the FIR. These include the Poisson noise caused by discrete 
unresolved sources and the clustering of galaxies, which to- 
gether represent the CIB anisotropics, the Galactic cirrus, and 
the instrument noise. Assuming that the noise is not correlated 
to the signal 



P(k) = y(k) [P sources (k) + Pdusik) + P c irrus(k)] + N(k), 



(2) 



Pdusik), and P C irms(k) are, respectively, the power spectrum 
of unresolved sources, the clustering, and the Galactic dust 
emission. The instrumental noise is represented by N(k), and 
y(k) is the power spectrum of the PSF of the instrument. 
The noise power spectrum, N(k) was computed using two 
independent maps of the ELAIS Nl field. At 100 //m, we used 
the different HCONs (see Sect. [24]). At 160 /mi, we used the 
even and odd BCDs to build two independent maps. The power 
spectrum of the difference in the two maps gives an estimate 
of N(k). To take account of the inhomogeneous cover a ge, w e 
applied the method of Miville-Deschenes & Laga che] p005| ). 
We subtracted the estimated N(k) from the raw power spectrum 
P(k). The level of the Poisson noise at large k is obtained by 
dividing the power spectrum by that of the PSF. 
The result is sho wn in Fig. [TT We found that P sources - 
9013 ± 100 Jy 2 /sr. |Lagache et aLH2007t found a slightly higher 
value of 9848+120 Jy z /sr with sources removed with a higher 
flux (200 mJy). In this study, sources are removed to a lower 
flux cut, which leads to a lower shot noise level. 
Statistical errors in power spectra are computed using mock 
signal plus noise maps that we analyzed with the same pipeline 
as the data. We derivedthe covariance matrix of this set of power 
spectra. Its diagonal terms give the errors in each P(k). Errors 
in the subtraction of the cirrus component are not statistical but 
systematic and due to errors in the emissivities dai (see Sect 
5.1 ). Errors in the power spectrum of the cirrus- subtracted map 
due to the spatial removal are on the order of (da) 1 x P(k), which 
is negligible compared to the statistical error. 
We have seen in Sect. 12.11 that there are uncertainties in the 
determination of the PSF. We took these into account by using 
several PSFs (ours, STiny Tim, and one computed from another 
Spitzer field, the CDFS) in the same pipeline. We derive again 
the covariance matrix of the set of power spectra. We add these 
errors in quadrature to statistical errors. 



where cr k = 0.065 ± 0.005 arcmin -1 at 100 jam, which corre- 
sponds in real space to a Gaussian function with cr = 1.8 ± 0.1 
arcmin. Sources are removed down to a 1 Oct threshold, follow- 



ing the algorithm described in Miville-Deschenes & Lagache 
2005 ). The source- subtracted map at 100 [im is shown in Fig. 



ma 

The FWHM of our data are given in Table 2.4 We convolved 
the maps when necessary in order to compare consistent data. 



4. The Galactic component 

In previous works, the Galactic contribution to the power spec- 
trum was either assumed to be proportional to k~ 3 (Lagache & 
Puget 2000 ), estimated usi ng IRIS at 100 urn o n large angular 
scales in the "&" 3 regime" (Lagache et al.||2007 ), or consi dered 



to be negligible (|Viero et al.|2009| ). |Matsuura et al.| ( |2011| ) used 
the |Schlegel et al. ( 1998| ) map at 100 fim as a cirrus template. 
This step was incorrect as, in the extragalactic fields, this map is 
highly contaminated (even dominated on some spatial scales) by 
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Fig. 11. The black line shows the raw power spectrum, P(k), 
of the map at 160 jum. The red line is the noise- subtracted 
power spectrum P(k) - N(k). The blue one represents the noise- 
subtracted power spectrum divided by that of the PSF (P(K) - 
N(k))/y(k). The dashed line shows the level of Poisson noise, 
which is 9 x 10 3 Jy 2 /sr. For clarity, the error bars are not shown. 



the anisotropies of the CIB, as demonstrated in Sect. 5^ 
The goal here is to examine the accuracy of the cirr us removal 
using Hi data (as done by Planck Collaboration et al. (2011b )), 
with the ultimate aim of obtaining CIB anisotropy maps on fields 
larger than ELAIS Nl, and even for a very large fraction of the 
sky with Planck data. In this section, we present our power spec- 
trum of the Galactic component using only 100 jim data on large 
scales. We then characterize the dust properties in our field. This 
is used in Sect. 5 to qualify the accuracy of the cirrus removal 
using Hi data. 



4.1. Constraints on the cirrus contribution using IRIS 100 \im 
data 

The power spectrum of the Galactic component i s first assumed 
to be a power law, as in Miville- Deschen es et al. ( 2007| ) 



s(k) 



(3) 



where Pq is the normalization of the power spectrum at 
ko = 0.01 arcmin -1 . The cirrus component dominates the power 
spectrum on large scales (k < 0.01 arcmin -1 ). The ten square 
degree field selected at 160 jim does not allow us to probe this 
regime reliably, therefore a larger map is needed. After remov- 
ing the sources, we fit the power spectrum of a 225 deg 2 IRIS 
100 fim map centered on ELAIS Nl deriving (3 = -2.53 ± 0.03 
and P = (4.93 ± .20) x 10 6 Jy 2 /sr. To compare this value to 
previous work (e.g. Lagac heet al.| ( |2007] ) in the Lockman Hole 
field), we need to rescale Po by the ratio of the cirrus bright- 
ness. We note that Po is indeed proportional to the square of the 
surface brightness of the cirrus B 2 ci = (#100 _ B C ib) 2 - If we 
take Bcib = 0.78 MJy/sr according to Laga che et al.| ( |2000| ), we 
1.25 MJy/sr. Using B La ^ che = 0.51 MJy/sr, we 



obtain B Us 

cirrus 

then have 



pLagache 
0,normalised 



pLagache 

r x 



B 



Us 



4.20 + 0.93 x 10 6 Jy 2 /sr(4) 



2 We therefore we did not compare their determination of the CIB 
with ours as their determination is highly biased owing to an incorrect 
cirrus removal. 



which agree with our measured Pq. 

The 100 yum power spectrum needs to be scaled by the ratio 
(#i6o/#ioo) 2 m order to be compared directly with the 160 
fim power spectrum ( |Roy et aL]|2010| ). This ratio depends on 
the cirrus physical properties and has thus to be determined 
sepcifically for our field. Unfortunately, this ratio is also scale- 
dependent: on large spatial scales, the brightness of the map is 
dominated by the cirrus, whereas on small spatial scales, CIB 
anisotropies dominate. We therefore differentiate the spatial 
scales larger than 95 arcmin from those smaller than 95 arcmin 
using a wavelet decomposition^] Moreover, we remove scales 
smaller than 6 arcmin as they are dominated by instrumental 
noise. We then perform a linear regression between the maps 
at 100 and 160 jjm for both spatial scales, yielding the ratio 
#ioo/#i60- To increase the statistics, we do a similar analysis 
in the Lock man Hole (Lockman et al. 1986| ), using the data 
published by Laga che et al.| (2007). Our results are listed in 
Table [2] They are compared with previous la rge-scale cirrus 
#ioo/#i60 colo r measurement s (Boulanger et al.|1996| ) and CIB 
colors from the |Lagache et al. (2003 ) empirical model of galaxy 
evolution. On scales <95 arcmin, our results are consistent with 



the CIB prediction of Lagac he et al.| ( |2003| ). However, there is 
a discrepancy on scales > 95 arcmin between our results and 
those of Boulang er et al.| ( |1996| ), which we attribute to real 
changes in cirrus properties from one field to another. Flagey 
[eTaLl ( [2009] ) found that B m /B m = 0.25 ± 0.01 across the 
whole Taurus complex, whereas l ocally, they foun d that the 
same ratio varies from 0.27 to 0.5. |Bot et al.| ( |2009| ) computed 
colors in several small regions of the sky surrounding nearby 
galaxies (-0.1 deg), and they also found varying colors from 
one field to another: from #i o/#i60 = 0-36 to #ioo/#i60=0-60. 
In both cases, they explained these differences by a variation 
of the interstellar radiation field and/or the abundance of very 
small and big grains. 

To compare the cirrus power spectrum obtained using IRIS data 
at 100 /dm with our 160 yum power spectrum, we rescale the 



100 yum power spectrum by B 



100,c 



/R 225deg* 
s /n lOO,cirn 



and then by 



(#i6o/#ioo) 2 = (1/0.35) 2 , the large-scale color that corresponds 
to the cirrus color in our field. We show this comparison 
in Fig. 12 and observe that they are in good agreement for 
k <0.02 arcmin -1 (i.e. in the "cirrus regime"). Unfortunately, 
the statistics of the large scales at 160 Lim are quite poor and do 
not allow us to quantify the quality of the argument. 



4.2. Dust colors and dust-Bi emissivities 

The dust that is heated by the interstellar radiation field and 
emits in the IR is mixed with neutral hydrogen. Thus, the IR 
emissi on of the cirrus is strongly correlated with the Hi 21 cm 
line. Boulanger & Peraulf| p988| ) showed that this correlation is 
tight at high Galactic latitudes at 60 and 100 jim. This correla- 
tion has o ften bee n used to study dust properties, for instance by 
|Boulanger et al. (1996]) who derived the dust spectrum associ- 
ated with Hi gas. 

In this section, we use GBT data at 21 cm to derive the FIR 
emission of the cirrus, which is then removed from our data at 
100 yum and 160 jim in Sect. [5] 

In the ELAIS Nl field, there are three distinguishable Hi veloc- 
ity components : the local component, an intermediate velocity 
cloud (IVC), and a high velocity cloud (HVC). We first compute 



IDL atrou algorithm 
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Field 



Bwo/Bu 



B\oo/B m >95' Bioo/B\eo < 95' 



Lockman Hole 
ELAIS Nl 



0.62+0.01 
0.50+0.01 



0.30+0.01 
0.35+0.01 



Cirrus color 




CIB colo r 
iLagache et al.|d2003b 



0.50 
0.27-0.5 
0.36-0.6 



0.64+0.01 
0.76+0.01 



0.65 



Table 2. Brightness ratios #i o/#i60 at scales larger and small er than 95 arcmin in the Lockman Ho l e and i n the ELAIS Nl field. 
The last four l i nes gi ve the cirrus and CIB colors according to |Boulanger et al. (1996), Flagey et al. ( 2009| ), Bot et al. ( 2009 ), and 
|Lagache et al. ( 2003| ), respectively. 




0.01 

k (arcmin -1 ) 

Fig. 12. Comparison of the cirrus power spectra computed with 
the IRIS 100 /mi map on large scales with that computed with 
the 160 fim map. Crosses : cirrus power spectrum at 100 yum 
computed on a 225 deg 2 100 yum map centered on ELAIS Nl 

rescaled by the ratio B™ 68 ? IB 22 ^ deg2 x (B m /B m ) 2 is shown 

J I00,cirrus' I00,cirrus v iDU/ iUU/ 

by crosses. Squares : power spectrum of ELAIS Nl computed at 
160 /mi. 



their integrated emission by adding all velocity channels with - 
14 km/s < V LS r < 43 km/s for the local, -79 km/s< V LSR <-14 
km/s for the IVC, and -163 km/s < V LSR <-79 km/s for the HVC. 
Second, assuming the optically thin case, we estimate their col- 
umn density using 



N H i(x,y) = 1.823 x 10 18 ^ T B (x, y, v)dv, 



(5) 



where Nhi is the Hi column density in unit of 10 20 atoms/cm 2 , 
Tb the brightness temperature, and v the velocity. 
The interstellar medium in the ELAIS Nl field seems to be dom- 
inated by neutral atomic hydrogen that reaches a peak Nhi ~ 
1.5 x 10 20 cm -2 in each of the three components (Fig. [8}. The 
brightness temperature of the Hi line is always < 8.9 K. Since 
molecular hydrogen, H2, begins to be see n in directions wher e 
Nhi > 2 x 10 20 cm" 2 and T b > 12 K (e.g.,|Gillmon et al.| ( |2006l ), 
Lisztj P002| ), |Lockman & Condon| ( [2005] )), it is unlikely that 
there are significant amounts of H2 in our field. Therefore, we 
can apply the decomposition following Miville-Deschenes et al. 
(2005) 



where lx is the infrared map, N l HI (x,y) is the column density of 
the /-the Hi component, a\ is the emissivity of component i at 
wavelength A, and Cx(x, y) is a residual term (offset + CIB). The 
correlation coefficients a\ are estimated using a x 2 minimiza- 
tiorj^] The error bars given by the IDL function are valid only if 
the noise of the lx maps is Gaussian and if the noise affecting 
Nhi is negligible. This may not be the case as the maps con- 
tain the IRIS or MIPS instrumental noise and CIB anisotropies. 
|Planck C ollaboration et al. (201 la) carried out Monte Carlo sim- 
ulations to estimate the errors in a\ for IRIS 100 and 60 yum and 
determined the coefficients by which they multiplied the error 
bars found assuming a Gaussian noise. We multiply our errors 
by these coefficients at 60 and 100 fim. For MIPS at 160 yum, 
we take the mean of the 100 and 350 jum coefficients determined 
by Plan ck Colla boration et al.| ( |2011a| ), as these coefficients vary 
only slightly with wavelength. They are on the order of eight. 
The emissivities a\ are computed at 60, 100, and 160 fim in 
the ELAIS Nl/MIPS field (i.e. ELAIS Nl field restricted to 
the MIPS coverage) and at 60 and 100 yum only in the ELAIS 
Nl/GBT field (i.e. the entire field covered by Hi data, see Fig. 
[8). Fig.[8]shows that ELAIS Nl/MIPS does not contain the IVC. 
Therefore, for this field, we use only two components, the local 
and the HVC, to avoid any additional noise. Our results are given 
in Table O 

The emissivities of the local component are in accordance with 
those found by Miville-Deschenes et al. ( 2005 ) at the three wave- 
length They used Hi observations from the GBT ([Lockman & 
[Condon 2005 ) to compute emissivities at 24, 60, 100, and 160 
fim in the Spitzer Extragalactic First Look Survey field. They 
have two IVCs whose emissivities are in agreement with ours at 
100 fim. However, there is a discrepanc y at 60 fim. For IVC1 and 
IVC2, |Miville-Deschenes et al.| ([2005} found B 60 /B m = 0.50 
and 0.34 respectivel y, whereas we have i?6o/ffi o o=0-30 . This 
value is in line with |Planck Collaboration et al!] ( [201 la| ), who 
found 0.23< B 60 /Biq <0.42 in 14 fields covering 800 square de- 
grees of the sky. They interpreted these differences as a contami- 
nation of the emission at 60 fim by non-equilibrium emission due 
to a higher relative abundance of very small grains compared to 
big grains at shorter wavelength. 

Moreover, the HVC in ELAIS Nl is detected at 160 fim. This 
confirms the detection of dust in some HVCs. Finally, we note 
that our results are also in very good agreement with those from 
Planck Collaborat ion et al. ( 2011a| ), who performed a similar 



analysis on the same field. 
Using emissivities, colors can be computed and compared to pre- 
vious _works_ 1 We find that #i6o/#ioo ~ 2.87 for the local compo- 
nent. [Botetal] ([2009]) removed the mean value of the CIB, 0.78 



Y J ^ i Hi(^y) + C A (x,y), 



(6) 



we use the IDL function regress 
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Fig. 13. The left panel shows the IRIS 100 fim map projected on the GBT map and convolved with the GBT beam. The emission of 
the cirrus (local+IVC+HVC) at 100 yum (the emissivity xN H i) is shown in the middle panel and the right one represents the IRIS 
100 /mi map convolved by the beam of the GBT from which we have subtracted the cirrus emission. 





Dust/Hi correlation coef in ELAIS Nl/MIPS 


Dust/Hi correlation coef in ELAIS Nl/GBT 


A (fim) 
a\ ocal this paper 
o^y c this paper 
a x this paper 


60 100 160 
0.175 + 0.056 0.87 + 0.14 2.46 + 0.43 

-0.004 + 0.015 0.034 + 0.031 0.31 + 0.08 


60 100 
0.175 + 0.014 0.877 + 0.028 
0.207 + 0.015 0.699 + 0.038 
-0.001 + 0.007 0.010 + 0.023 


af al MAMD 2005 
a 1 / 01 MAMD 2005 
a 1 / 02 MAMD 2005 
a H x vc MAMD 2005 


0.16 + 0.02 0.80 + 0.08 1.7 + 0.02 
0.35 +0.04 0.70 +0.09 2.7 +0.4 
0.31+0.04 0.9+0.1 1.4+0.4 
0.05+0.01 0.055+0.015 0.8+0.1 




af aL Planck 2011 
a[ vc Planck 2011 
a^ vc Planck 2011 




0.166+0.011 0.862+0.033 
0.213+0.012 0.723+0.036 
-0.001+0.007 -0.009+0.022 


< residue > (MJy/sr) this paper 
CIB mean (MJy/sr) Bethermin 201 1 


0.58 0.76 a 
0.30+0.01 0.63+0.01 





Table 3. Emissivities and CIB levels. Columns 4, 5, 6 : Emissivities (in units of MJy/sr (10 20 H atoms) -1 cm 2 ) computed in ELAIS 
Nl/MIPS at 60, 100 and 160 jim compared to those of the literature. Columns 7, 8 : Emissivities computed in Nl/GBT at 60, 100 
fi m. Uncertainties are 1 cr uncertai nties that take into account the statistical variance and the instrumental noise . MAMD 2005 stands 
for[Mi ville-Deschen es et aL| ( |2005| ). Their emissivities are for the Spitzer XFLS field. Planck 2011 stands for Planck Collabo ration! 
|et al.| ( |20lTa ). Their emissivities are for the ELAIS Nl/GBT field. The second last line gives the mean values of the residual maps. 
At 160 fim, the mean value has been corrected for the offset determined in Sect. 6.1, and gives the level of the CIB (although see 
Sect. 6 for a more accurate determination). At 100 yum, the mean value is strongly contaminated b y residual zodical emiss ion (see 
Dole et aL] ( |2006| )) thus does not give the CIB level. The last line gives the CIB coming fr om th e Bethermin et al. ( 2011| ) model. 
{a) This value has been corrected from an offset that is present in the scanning map (see Sect. 6.1 ). 



MJy/sr from their < #ioo >, and they observed an increase in 
the ratio #i6 /#ioo with < #i o > meaning that brighter regions 
are colder. We subtract the same CIB level of 0.78 MJy/sr from 
our < #100 > , even though we get a mean CIB residual of 0.58 
MJy/sr (see Tab[3]). This results in < # 10 o >= 0.4 MJy/sr, which 
is much lower than any of their < #ioo >. They kept only re- 
gions with < #100 »2.5 MJy/sr, in order to be dominated by 
variations in cirrus emission (and not be contaminated by CIB 
anisotropies). Their fields have a typical #i6o/#ioo of 2, which 
is very close to our value. Our much lower < #ioo > implies 
that #i6o/#ioo may reach a plateau at a value of ~ 2 MJy/sr for 
< #ioo > lower than 2.5 MJy/sr. The value of #i6 /#ioo can also 
be compared to those given in Table [2] which were computed 
with a linear regression between the two maps at large scales. We 
achieve an excellent agreement, finding #i o/#i60 = 0.35+0.01 
in the two cases. 

The CIB should be the only astrophysical component contained 
in the residual map. To test this, we can check whether the resid- 
ual mean value agrees with the mean value of the CIB. This is 



the case at 160 yum, after the map is corrected by an offset term 
(see Sect. 6.1 ). At 100 yum, there is a discrepancy between the 
prediction and our result, which we attribute to residuals of zo- 
diacal light in the map. The IRIS data were indeed calibrated in 
the DIRBE data, which contain a residual zodiacal emission that 
leads to an overestimate of the CIB (see Dol e et al.| ([2006 )). At 
60 yum, we do not try to measure the mean value of the resid- 
ual because the CIB level is on the order of the residual of the 
zodiacal light at this wavelength (Renault et al. 2001 ). 



4.3. Dust temperatures 

The IR/Hi emissivities give constraints on the dust temperature. 
We assume that the emission of big grains at thermal equilibrium 
with a radiation field is a modified black body 



Iy = TyBy(T BG ), 



(7) 



where B v is the Planck function, T BG the big grains' equilib- 
rium temperature, and r y the optical depth. It can be expressed 
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as r y = Nhi€ v with e v the dust emissivity per H atom that is 
usually assumed to be a power law e v = eo(v/vof where ft is 
the emissivity spectral index. Following Boulan ger et al.| ( [T996] ), 
we assume that J3 = 2. At wavelengths where the emission is 
dominated by big grains, the IR/Hi correlation coefficients can 
be written as a v = e v B v (T BG ). We use the infrared-Hi correlation 
coefficients at 100 and 160 jim to estimate the big grains temper- 
ature T l BG of each Hi component assuming that the contribution 
of very small grains in the far-infrared is negligible. 



1 o 7 



'100 



£ y (100//m, 7^) /160 



(8) 



Bot et al. 



(2009), 



a[ 60 B v (160fim,P BG )\100 

We find T 1 "™ 1 = 15.9 ± 0.2 K. This agrees with 
who found dust temperatures in the diffuse medium ol between 
15.7 K and 18.9 K and with |Miville-Deschenes e7alj ( |2005] ) who 
found T l °™ 1 between 16.3 and 18.8K. We do not detect the HVC 
at 100 yum, hence we cannot determine its temperature. However, 
we compute a limit to its temperature assuming a 3<x limit for the 
detection at 100 fim. We get T BG C < 15K, which agrees with 



BG 

Miville-Deschenes et al. (2005k who found a value of T BG C be- 
tween y.y K and ll.b K. The HVC is colder than the local dif- 
fuse medium. This is consistent with a lower radiation field than 
in the solar neighbourhood owing to the distance of the HVC. 



5. Removal of the Galactic component 

5.1. Removal of the cirrus either spatially or by considering 
the power spectra and error bars 

The cirrus contribution removal can be done in two ways: spa- 
tially in the maps or by subtracting the cirrus power spectrum 
from that of the infrared map. We first need to know whether 
these two methods are equivalent. We consider the power spec- 
trum of the CIB, PciB(k), as an estimator for this test, and com- 
pute its variance for the two methods in a naive approximation 
of Gaussianity. We have two maps A and B, containing, respec- 
tively, CIB + Galactic component (GC) and only the GC. In 
Fourier space, we have 4 = a^ciB + &k,GC an d 4 = cl^gc? 
where are the Fourier coefficients. The power spectrum can 
be written 



PcibQc) 



Z< a k,CIB > 



(9) 



The mean values of the power spectrum of the CIB anisotropies 
are equal for both methods but not their variances. For the spatial 
removal, 



Var(PciB) = < (4 - a B k f > - < (4 - a*) > 2 
= Var(PciB), 



(10) 
(11) 



whereas for the removal at the power spectrum level a correla- 
tion term between the cirrus power spectrum and that of the CIB 
appears, in addition in the variance of the CIB power spectrum 



var(PdB) = < (44 - 44f >-< (44 - 44) > 

= Var(PciB) + Var(PciB) x (4?. 



(12) 

(13) 
(14) 



Thus, from this simple argument we expect error bars to be 
smaller after the spatial subtraction. We check this by perform- 
ing simulations using mock data. We generate signal plus noise 
maps and analyze them in both ways. Error bars are indeed 
smaller in the spatial subtraction case. The template (spatial) 
subtraction removes each moment of the statistics, whereas the 
second method only removes the moment of the power spectrum. 



1 0° r 



1CT r- 




0.01 



k (arcmin 1 ) 



Fig. 14. The black continuous line with the crosses and bold er- 
ror bars show the power spectrum of IRIS 100 jim and the red 
dotted line with squares and thin error bars represents the power 
spectrum of the emission of the cirrus at 100 /mi calculated from 
the Hi map. The discrepancy observed for k>0.02 arcmin -1 is 
due to the clustering of SB galaxies (see Sect.|6.2|). 



With cirrus 
Cirrus free 




Fig. 15. Red crosses show the CIB anisotropies power spectrum 
computed at 100 /mi. The c ontribution of the cirrus has been re- 
moved following Sect. 5.2 The power spectrum of the noise has 
been subtracted and the resulting power spectrum has been di- 
vided by the power spectrum of the PSF. The black line is the 
power spectrum of the map that contains the cirrus contamina- 
tion. The dashed line shows the level of the Poisson noise which 
is 8.7xl0 3 Jy 2 /sr. 



5.2. Contribution of the Galactic component to the power 
spectrum at 100 j^m 

We start by comparing the cirrus power spectrum obtained at 
100 jim to that computed with Hi data in the ELAIS Nl/GBT 
field. To construct a map of the cirrus emission at 100 /im, we 
take the sum of each component colu mn density weighted by the 
emissivities determined in Sect. 14.21 : 



d _ Aocal sy A jlocal , ^IVC v \7 IVC _i_ ™ HVC n/ A7 hvc 
#100 - ^100 XJ Vmt +^mn X JN „ T + Of im X 7V ft 



"100 



'100 



(15) 



The left panel of Fig. 13 shows the map at 100 jim projected 
onto the GBT map and convolved with the GBT beam, the 
middle panel shows #100 computed from 21cm data. The right 
panel shows the 100 jim cirrus-free residual. We can clearly 
see that most of the contribution of the cirrus is removed and 
that the residual map contains only the CIB. This is even 
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more visible in Fig. 14 which compares the 100 fim power 
spectrum to that of #100 computed from 21 cm data. They are 
in very good agreement for k <0.02 arcmin -1 . The discrepancy 
observed on smaller scales (k>0.02 arcmin -1 ) is exp ected, as the 
clustering of SB galaxies dominates (see Sect |6.2|). Therefore, 



this shows that using of IR/Hi emissivities is a reliable method 
for measuring the contribution of the Galactic component in the 
far-infrared. 

We can now subtract this cirrus emission from the origi- 
nal map to get the CIB anisotropy map. Since the FWHM of 
the GBT is much larger than that of IRIS, the 100 yum ELAIS 
Nl/GBT map must be convolved with the GBT beam (9.1 
arcmin). Since the cirrus contamination is primarily on large 
scales, we can create a hybrid power spectrum using the power 
spectrum of the CIB anisotropies obtained by removing the Hi on 
large scales, and that of the original map on small scales. We use 
k = 0.05 arcmin -1 for the transition between the two power spec- 
tra. We note that the power spectrum of the CIB anisotropies on 
large scales must be divided by the power spectrum of a PSF 

with a FWHM = ^(FWHM 2 GBT - FWHM 2 IRIS ) to take into ac- 
count the convolution by a larger beam. The noise map is com- 
puted following Sect. [3] and the PSF is that of Sect. [2] We sub- 
tract the noise power spectrum from the raw one and divide it 
by the power spectrum of the PSF. We present the power spec- 
trum as well as the cirrus-free one on Fig. [15] We clearly see the 
amount of power due to the cirrus that has been removed. We 
derive a Poisson noise of 8690± 190 Jy 2 /sr, which disagrees with 
Miville-Deschenes et al] ( |2002| ) who found ~5.8xl0 3 Jy 2 /sr, al- 
though the bright point sources are removed in the same way. 
The main difference between both studies is the removal of the 
cirrus contribution. They fitted the large-scale part of the 100 jim 
power spectrum and removed the fit on all scales. When we carry 
out the same analysis, we find a Poisson noise ~ 5xl0 3 Jy 2 /sr, 
in agreement with them. We conclude that they overestimated 
the cirrus contribution by removing a power law on all scales 
applying a fit only at large scales. The power-law fit is indeed 
contaminated by CIB anisotropies. 

5.3. Contribution of the Galactic component to the power 
spectrum at 160 jim 



Using the same method as in Sect. 5.2 we remove the cirrus 
emission from the ELAIS Nl/MIPS map at 160 jim. We com- 
pute the hybrid power spectrum with a cut at k = 0.05 arcmin -1 . 
Fig. [16] shows the total power spectrum (black) and the cirrus- 
free one (red). We clearly see the difference only on the largest 
scales available with this map. We also plot our fi t to the 
clustering power spectrum from Laga che et al.| p007| ) in blue, 
and the shot noise level in black. The green line shows the 
sum of the two-component fits. We also compare our resulting 



power spectrum to that of Lag ache et al.| ( |2007| ) in Fig. [T7 



These power spectra are in very good agreement for k > 0.03 
arcmin -1 , where they are dominated by the CIB anisotropies 
(both clustering and Poisson noise). On scales < 0.01 arcmin -1 , 
there is more power in the Lagach eet al.| ( |2007] ) power spectrum 
because it contains the cirrus contribution (the blue dashed line 
is their estimate of the power spectrum of the cirrus). We can 
see that using Hi data, we are able to extend the measurement 
of the correlated fluctuations to large scales. This shows that 
making use of Hi data at 21 cm is an efficient way of removing 
the contamination of the Galactic component. 



Clustering + Poissonian fit 
Poissonian fit 
Clustering fit 

With cirrus 
Cirrus free 



0.01 0.10 
k (arcmin -1 ) 

Fig. 16. MIPS power spectrum at 160 yum. The dashed black line 
shows the power spectrum of the raw map (noise subtracted and 
divided by the PSF power spectrum). The red line represents the 
power spectrum of the map to which the cirrus emission has been 
subtracted. The blue dash-dotted line shows our fit to the power 
spectrum of the clustering and the black horizontal dotted line is 
the fit to the Poisson noise level. The green line shows the sum of 
the clustering and Poisson components. Error bars are not shown 
for display purposes. 
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Fig. 17. Black squares are the resulting power spectrum of the 
CIB anisotropies computed by removing the cirrus contamina- 
tion using Hi data in the ELAIS Nl field. Red crosses repre- 
sent the Lagac he~et al.| ([2007) power spectrum computed in the 
Lockman-Hole field. It includes both CIB and cirrus anisotropies 
and the blue dashed line shows their estimate of the cirrus power 
spectrum. 



6. Measurement of the CIB mean levels at 100 /zm 
and 160 jim 

6.1. CIB mean estimate with Total Power Modes at 160 \im 

We can also use our Hi data and emissivity measurements to 
compute the absolute level of the CIB at 160 pm. We consider 
the two TPMs of Spitzer archival observations (26961920 & 
26962432) of the ELAIS Nl field that were designed to cross- 
check the calibration of the diffuse emission at 160 pm. Even 
with a cryogenic telescope such as Spitzer, there is a small com- 
ponent of thermal emission at longer wavelengths that contam- 
inates the standard photometric observations. The TPM mode 
by-passes the effects of this spurious radiation by comparing the 
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Fig. 18. Positions of the total power mode observations on the GBT map with their Hi velocity profiles. 



emission of the target (sky) with that of an internal dark to pro- 
vide an absolute measurement (see MIPS Handbook § 3.1.12). 
This mode was designed precisely to observe relatively faint ex- 
tended emission r egion s. The TPM observations that we used are 
discussed in Sect. 12.21 

There is an HVC in TPM 1 and only the local component in TPM 
2 as shown on Fig. [18] We first compare the MIPS scan map val- 
ues to the TPMs. By first subtracting the map values from those 
of the TPMs (see Table [4|, we determine the offset of the scan 
map from each TPM position. They are in good agreement in the 
two regions with an average offset of -2.05+0.24 MJy/sr, which 
has no consequences for the power spectrum estimate as well as 
on the CIB level determination that follows. 
Making use of the TPM values and the emissivities previously 
computed, we determine the absolute level of the CIB at 160 



/mi. After subtracting the zodiacal light, the components of the 
TPM are 



TPM - Zodiacal = Cirrus + CIB, 
where the cirrus brightness is 



Aocal ■ 



jlocal 



JVC . 



rIVC 



HVC 



tHVC 



(16) 



(17) 



and a are the values computed across the whole ELAIS 
Nl/MIPS field, as listed in Table [3] We use emissivities calcu- 
lated over the entire field rather than over the TPM regions alone 
since they are far too small to obtain an accurate measurement. 
The zodiacal light, which has a value of B zo< n - 0.83+0.12 
MJy/sr, was previously subtracted. It wa s estimated using th e 
Spitzer background model (Reach 2000; Reach et al. (1995); 
Kels all et al. ( |1998| )). With our emissivities, we can compute 



the cirrus brightness of both TPMs (following Eq. 17) and thus 
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compute the CIB level (following Eq. 16). Our results are listed 
in Table [4] The last column gives the CIB levels: they are in 
excellent agreement, even though they were computed in two 
regions with completely different cirrus contaminations. We get 
an average value of B m = 0.77+0.04+0.12 MJy/sr for the CIB 
level at 160 yum. The first error is statistical, and the second 
is systematic and dominated by the error in the zodiacal light 
removal. We do not include the calibration uncertainties (that 
are on the order of 12%, ([Stansberr y et al.|2007| )) in our quoted 
errors. The errors are dominated by those of the zodiacal light 
model. 

We can compare the value of the CIB level determined 
with the TPMs to the mean of the residual map we obtained 
in Sect. |4.2| as the only astrophysical component that should 
be present in this map is the CIB. We have first to correct the 
map for the offset. We use 2.05 MJy/sr and find B resi d U ai = 0-76 
MJy/sr. This value is in very good agreement with that ob- 
tained using the TPMs. They are very close to the predicted 
CIB le vel of B 160 = 0.63 + 0.02 MJy/sr by |Bethermin eTaL 
( 2011) and to the last Spitzer determi nation using very deep 
number counts: pethermin et al.| ( |2010| ) found B^q = 0.78^ ^ 
MJy/sr. Moreover, our value is also in very good agreement with 
that of |Berta et al.| ( [20TT] ), who carried out a P(D) analysis in 
addition to a stacking one (of sources detected at 24 fim) in 
Herschel/PACS data at 160 fim to derive the differential number 
counts. Extrapolating the counts down to very faint fluxes using 
a power law, they obtained B^q = 0.72^-^ MJy /sr. We can also 



compare our determinations with Juvela et al. (2009), who de- 
rived the CIB using ISOPHOT data. They computed linear fits 
between FIR and Hi data, by disregarding several velocity com- 
ponents in Hi data and for much smaller fields than ELAIS Nl 
(roughly 25 times smaller). They found B m = 1.08 + 0.32+0.30 
MJy/sr in the range 150-180 yum, where the first error is statis- 
tical and the second systematic. Although compatible within the 
error bars, our determination points to a lower value of the CIB. 
Our CIB determination benefits from having Spitzer and GBT 
data that cover a large field, allowing for a more robust measure- 
ment of the cirrus contamination and thus of the CIB. 
We can now combine our CIB mean level measurement at 
160 yum with our CIB anisotropy measurements to compute the 
CIB at 100 fim. 

6.2. CIB mean estimate at 100 fim 

We combine the cirrus-free power spectra measurements at both 
100 yum and 160 yum with the CIB mean at 160 yum to derive the 
CIB mean level at 100 fim, following 



(T 



CIB 
100 fim 

n \00fim 



CIB 
u 160 fim 

t>CIB ' 
U 1 60 fim 



(18) 



where cr is the rms fluctuation in the CIB that can be computed 



-CIB 



using the measured CIB power spectrum. Knowing cr\ Q0 m , 



r CIB 

160yura' 



160 fim' 



we can compute B^^ 



following Eq. 18 



We checked using the model of Penin et al. (2012 ) that Eq.[TT 
valid. 

The two CIB power spectra determined previously at 100 yum 
and 160 jjm were measured with a different flux cut for the bright 
sources removal. Thus, we first recompute the power spectra for 
the two maps without masking any bright point sources, in or- 
der to have the same point source contamination (i.e. consis- 
tent Poisson noise levels). We then perform the same analy- 
sis as previously (removal of instrumental noise and cirrus, and 



division by the power spectrum of the PSF). We compute the 
100/160 color of the clustering power spectra using the largest 
common scales available (0.006< k <0.2 arcmin -1 ), obtaining 



r CIB 

lOOfim 1 ^ \60fim 



lo\£ = 0.3 1 + 0. 1 . Multiplying the CIB at 160 yum de- 



rived in Sect. 



R CIB 

n ioo 



6.1 



by of 7 * lcr CI] ? leads to a CIB at 100 fim of 



100 fim 1 160 fim 

0.24+0.08+0.04 MJy/sr, where the first error is statis- 
tical and the second one is systematic. We note that this value 
is lower than the mean value of the residual map obtained af- 
ter the removal of the cirrus component (see Table [3]), owing 
to residual zodiacal emission in the IRIS map. An empirical 
correction o f the |Lagache et aL] ( |2000| ) DIRBE measurements, 
provided in |Dole et al.| ( |2006| ) gives 0.48+0.21 MJy/sr. Our 
determination points towards a lower CIB value. By stacking 
Herschel/PACS maps at 100 fim at the positions of all 24 fim 
galaxies (S(24)>20yuJy), |Berta etaL] ( |2010| ) measure d a CIB sur- 
face brightness of 0.25+0.02 MJy/sr. Subsequently, Berta et al. 
( 201 1 ) improved on these measurements by carrying out a P(D) 
analysis that allowed them to reach lower fluxes. Extrapolating 
their differential number counts using a power law, they obtained 
0.42 + ^ MJy/sr, which is barely compatible with our 



pCIB _ 
100 vy '" r ^-0.06 

determination. 



7. Conclusion 

We have presented a new method for removing Galactic cirrus 
contamination from the power spectrum of CIB anisotropies by 
using an independent tracer of this cirrus, the Hi 21 cm data. We 
have computed the far-IR emission of each velocity component 
of the cirrus and removed it spatially from the maps. The resid- 
ual map is thus a map of the CIB anisotropies. We have applied 
this method to MIPS dat a at 160 /im in t he EL AIS Nl field and 
recovered the results of Lagache et al. ( 2007| ) on intermediate 
spatial scales where starburst galaxy clustering intervenes. They 
had derived a linear bias b ~ 2.4 with MIPS data at 160 yum, 
which probes mainly galaxies around z ~ 1. We similarly ap- 
plied this method at 100 yum and detected for the first time the 
correlated anisotropies at this wavelength. An analysis of these 
CIB anisotropy power spe ctra was not an obje ctive of this paper, 
and we refer the reader to Penin et al. ( |2012| ), who presented a 
clustering model of star- forming galaxies. 
We have shown that the cirrus removal using Hi data is the most 
efficient available method for removing the cirrus contamina- 
tion accurately. It has also been successfully applied to Planck 
data ( [Planck Co llaborati on et al.||2011b| ) and will be applied to 
Herschel data. We caution the reader th at it is incorrect to use 
the 100 fim map of Schle gel et al.| ( |1998] ) as a cirrus tracer as this 
map contains the anisotropies of the CIB. 
We have used absolute measurements of the brightness on small 
regions of the sky (total power modes) with different cirrus con- 
tributions to derive the CIB level. Making use of our emissiv- 
ities, we were able to compute the cirrus brightness in two of 
these regions, and derive the CIB level at 160 yum. We found 
B m = 0.77+0.04+0.12 MJy/sr (the first error is statistical and 
the second systematic). In addition, using our measured CIB cor- 
related anisotropies at 100 yum we computed the CIB anisotropy 
color, #ioo/#i60 = 0-31 +0.1. This color measurement is free 
of cirrus and zodiacal light uncertainties, since the former had 
been removed from both power spectra and the latter is a con- 
stant that has no influence on power spectra. Using this color 
and the CIB measured at 160 fim, we derived the CIB at 100 fim 
#ioo = 0.24 + 0.08 + 0.04 MJy/sr where the first error is statisti- 
cal and the second systematic. These CIB measurements are the 
most accurate measurements based on a diffuse emission anal- 
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Name 


Coordinates 


TPM value 


N Hi, local X a local 


Nm,Hvc x &hvc 


CIB level 


map 


offset 


1 


RA=245.20 Dec=54.98 


2.59+0.12 


1.40+0.03 


0.43+0.01 


0.75+0.04+0.12 


4.78+0.12 


2.20+0.24 


2 


RA=243.28 Dec=55.13 


2.23+0.12 


1.43+0.03 


0.01+0.01 


0.78+0.04+0.12 


4.14+0.12 


1.91+0.24 


average 










0.77+0.04+0.12 




2.05+0.24 



Table 4. TPM observations at 160 jim. The second column gives the coordinates of the TPMs, the third gives the brightness of the 
TPMs (zodiacal-light subtracted). Infrared emission of the local and of the HVC are in columns 4 and 5. The sixth column gives 
CIB levels at 160 jim. The first error is statistical and the second is systematic. It is due to the error on the estimation of the 
zodiacal light. The offsets between the scan map and the TPMs are given in the last column. All values are in MJy/sr. 



ysis. Comparing those measurements with an extrapolation of 
galaxy number counts, we have found no evidence of an un- 
known contribution to the CIB, in contrast to |Matsuura et al.| 
fl20TTT ). 
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